Temporal behaviour of global perturbations in compressible axisymmetric 

flows with free boundaries. 



V.V. Zhuravlev 1 , N.I. Shakura 1 ' 2 

*1) Sternberg Astronomical Institute, Universitetskii pr. 13, Moscow, 119992 Russia 
2) M ax-Planck- Institut fur Astrophysik, Karl-Schwarzshtld-Str. 1, Postfach 1311, D- 8 5741 G arching, Germany 



key words: accretion, accretion disks - hydrodynamics - instabilities - turbulence. 



Abstract 



o 
O 
in 



Oh 
i 



00 1 

o : 
o . 

The dynamics of small global perturbations in the form of linear combination of a finite number of non- 
axisymmetric eigenmodes is studied in two-dimensional approximation. The background flow is assumed to be 
an axisymmetric perfect fluid with the adiabatic index 7 = 5/3 rotating with power law angular velocity distribu- 
tion Q oc r~ q , 1.5 < q < 2.0, confined by free boundaries in the radial direction. The substantial transient growth of 
acoustic energy of optimized perturbations is discovered. An optimal energy growth G is calculated numerically for 
a variety of parameters. Its value depends essentially on the perturbation azimuthal wavenumber m and increases 
for higher values of m. The closer the rotation profile to the Keplerian law, the larger growth factors can be obtained 
but over a longer time. The highest acoustic energy increase found numerically is of order ~ 10 2 over ~ 6 typical 
Keplerian periods. Slow neutral eigenmodes with corotation radius beyond the outer boundary mostly contribute 
to the transient growth. The revealed linear temporal behaviour of perturbations may play an important role in 
angular momentum transfer in toroidal flows near compact relativistic objects. 

, c3L Introduction 

£\Hydrodynamical equations allow the existence of laminar axisymmetric toroidal flows with free boundaries in the 

^vicinity of a gravitating body. The fundamental question is whether the accretion is possible if small perturbations 

of some sort are present in the medium, i.e whether the perturbations will grow and affect the background to allow 

angular momentum transfer. Astrophysically, it is interesting to know if such a transition exists in flows with an 

^Lalmost Keplerian rotation. In the last case setting free boundaries leads to the hypersonic character of the basic flow 

with Mach number M » 1. 
GO 

, The stability of toroidal flows was studied extensively in the 1980s and 1990s by means of the traditional eigen- 
OQvalue analysis which has been widely used since the well-known investigations of Lord Rayleigh. The growing non- 
^Spdsymmetric definite frequency modes were discovered in the barotropic tori with a positive radial gradient of the 
^Specific angular momentum. This phenomenon was cold the Papaloizou-Pringle instability in honour of the authors 
^vho published their work in 1984. In numerous subsequent papers this sort of instability was studied in two- and 



hree-dimensional approximations by both analytical and numerical methods [25J [HI US El H3 ESI HH [HI EI] • It was 
^found that there are growing surface gravity and sound waves. The general picture of the instability is constructed 
■ oy these modes which grow due to the mutual coupling as well as by the Landau mechanism, i.e. due to interaction 
with the background flow 26, 28 . However, as the rotation profile approaches the Keplerian law the gravity surface 
waves become damping while the sonic instability ceases because of very small increments [HI EH] ■ This fact actually 
dropped the subject since the near-Keplerian rotation is of the most interest in connection with the problem of angular 
momentum transport in the accretion disk theory. Besides, the growing linear modes can exist owing to additional 
degrees of freedom in the flow. For instance, [7] and [36] claimed that stratorotational instability emerges in vertically 
stratified rotationally stable flows in the shearing sheet approximation. 

However, instability of axisymmetric flows can be considered from a quite different viewpoint. It is closely associated 
with the formulation of the initial value problem for perturbations. Indeed, to examine the stability thoroughly one 
must show that any initial perturbation will not ever affect the background flow. However, the absence of growing 
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eigenmodes guarantees only the asymptotic stability and tells us nothing about the behaviour of perturbations at 
finite time intervals. The generalized approach to the hydrodynamical instability was realized and declared in the 
most distinctive form among the fluid physicists in the middle 1990s [301 OH]- It is often called the nonmodal stability 
theory and concentrates on the initial value problem, partially on searching for optimal initial perturbations that can 
exhibit large transient growth even in asymptotically stable flows. Yet the role of this linear mechanism is usually 
emphasized in connection with the transition to turbulence |18j . According to the so called bypass concept the linear 
dynamics is responsible for extraction of the perturbation energy from the background. It is supposed that nonlinear 
processes only redistribute energy between the modes. The general description and some references to hydrodynamical 
literature on this subject can also be found in [35) and in application to astrophysical problems in the more recent 
paper [5]. 

Mathematically, the evolution of small perturbations is governed by a linear dynamical operator which in general 
can be normal or non-normal. In the first case eigenmodes are orthogonal and the cigenanalysis gives the complete 
picture of the temporal evolution of perturbations. Strictly speaking, the largest possible growth of perturbations is 
determined at any moment by the highest eigenvalue [3] . In the case of non- normal dynamical operator eigenmodes 
become non-orthogonal and the substantial temporal growth of perturbations is possible even if there are no growing 
eigenmodes [3"U] . 

The nonmodal approach to local perturbations was first pointed out also in astrophysical literature in [21j and 
|12j . The local framework allows us to perform the spatial Fourier transform and thus to study directly the temporal 
behaviour of perturbations. This direction was developed then in a number of papers [U, 5, 33, 38, [Tl |2"31 134) . 

Another branch of astrophysical investigations that use nonmodal approach is devoted to global dynamics in tur- 
bulent accretion disks. For instance, [5] considered the dynamics of small stochastic perturbations in an incompressible 
Keplerian disk and found the resulting steady-state coherent structures enhancing the angular momentum transfer 
outward. It is also worth mentioning the work |37j which studied the transient growth of axisymmetric perturbations 
in a thin accretion disk with account of compressibility and vertical structure of the flow. 

Here we concentrate on the global two-dimensional dynamics of the bounded laminar axisymmetric flows which 
were studied actively 15-20 years ago. We do not introduce viscosity since its action is negligible over the first ~ 10 — 50 
rotation periods but we include compressibility which is necessary to do since the sound speed is always comparable 
or less than the shear velocity difference when we are dealing with free boundaries. First we briefly describe the 
solution of the boundary problem, i.e. the calculation of eigenmodes. This will be followed by the discussion of the 
optimization algorithm which allows us to determine the specific group of eigenmodes with the strongest transient 
growth. This is actually the standard mathematical procedure that has been used by many authors. Then we present 
the results and make the conclusions. 

2 The Basic Flow and the Boundary Problem for Perturbations. 

We neglect here the vertical structure of the flow and consider an idealized barotropic cylindrical configuration rotating 
with the angular velocity power profile Q — Qo r ~ q (here fio is the Keplerian frequency at ro where the pressure 
attains its maximum). To make strict conclusions about the dynamics of real toroidal configurations in the vicinity 
of a gravitating body one must use full three-dimensional approach. In spite of that it is well known that for a true 
barotropic flow the problem can be solved using two-dimensional equations with a two-dimensional polytropic index. 
This is proved analytically for A >> h, where A is the wavelength of perturbations and h is the vertical scale of 
the basic flow (see [H] for details). Moreover, numerical studies have shown that the perturbed flow behaves almost 
like in two dimensions even for A ~ h 119] . The authors themselves explained that by the form of the Reynolds 
power which does not contain the vertical velocity perturbation in the case of true barotropic configurations. Since we 
are interested in the dynamics of an arbitrary linear combination of the flow eigenmodes with exponential temporal 
evolution cx exp(-iujt), we should start with solving the boundary problem. For the flow considered here this was 
done first by |13j and reexamined for growing perturbations in our previous paper [39] . The spectrum consists of 
an infinite number of acoustic modes and may be conventionally divided into three parts. The main part is the so 
called corotation interval, corresponding to the mode corotation radius (i.e. the point where u) = mil) inside the 
flow. The corotation interval is given by (mil (r 2)1 mf2(ri)), where r\ and ri are the inner and outer boundary points 
respectively. Owing to modes coupling and the Landau mechanism 28 , 32 this part of the spectrum contains damping 
and growing modes which have been of the main interest in previous investigations in the framework of the traditional 
eigenvalue approach. However, in flows with almost Keplerian rotation acoustic eigenmodes exhibit very slow growth 
thus weakly affecting the basic flow. The other two parts of the spectrum correspond to the modes with the corotation 
radius inside the inner boundary or beyond the outer boundary of the flow. As in this case the energy exchange 
with the background medium is impossible, they are pure neutral. In the present paper we restrict ourselves to study 
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the latter kind of eigen-solutions both with the growing eigenmodes from the corotation interval. This restriction is 
done because the Lin rule [20) forbids to calculate damping and neutral eigenmodes inside the corotation interval by 
integration along the real axis which has been used in our numerical method in the previous paper [39] . 

The eigen- or, alternatively, definite frequency modes have the standard form: 5h(r,tp,t) = h{r)e~ l ^ wt ~ mip ^ . Here 
to = 1,2,3, ... is the azimuthal wavenumber and u> — oj r + iu>i is generally a complex frequency with real part (also 
called the pattern speed) and imaginary part (also called the increment). Components of the perturbations velocity 
field (8v r , Sv v ) have the same form and it is easy to obtain the relations between v r , v v and h which we omit here but 
will be necessary below. 

The equation for eigenmodes in a barotropic flow can be easily written for h: 



where 
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— (£lr 2 ) is the epicyclic frequency 

r dr 



D = k 2 - (w - mtt) 2 

p, a 2 are the density and squared adiabatic sound speed profiles of the basic flow, respectively, which can be expressed 
in terms of the enthalpy h — a 2 /(7 — 1). In turn, h{r) has to be found from the basic flow stationarity expression: 



where il g = -J GM/r 3 is the Keplerian frequency. 
The boundary condition is: 
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where Ah, Ap are Lagrangian perturbations of the enthalpy and pressure. 
Expression © must be rewritten for h: 
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It turns out that ([3]) is equivalent to the regularity condition for the solution at the boundary points, which must 
be imposed since coefficients of eq. (TTJ become singular at r\ and r%- 

Eq. |T]) with b.c. © was solved for a variety of parameters in [5S]. As an example, the part of the spectrum 
for q = 1.58, w = 300%, 7 = 5/3 and to = 5 is shown in Fig. 1. In this figure the first several neutral modes with 
uj < mQ,{r-2) and some growing modes from the corotation interval with uji > 10~ 4 are displayed. Here and below 
all frequencies are given in units of Slo and time in units of f^ 1 ; we also introduce w = (rz — ri)/ra x 100%, the 
radial size of the flow in percents. We also use the Mach number defined as M = (f2(ri)ri — Q(r2)r2) / a(ro) though 
the presence of free boundaries where p = a = and non-trivial distribution of a(r) makes M to be not so obvious 
parameter to describe the compressibility of the flow. In real three-dimensional flows (i.e. tori) around a gravitating 
body M approximately equals to 1(ri — n)/h, which is the ratio of the torus radial size to its half- thickness. 



3 The Optimization Problem 

3.1 Class of perturbations and the growth measure choice 

Let us consider a linear subspace of the solutions of the perturbation equations, which is represented by the linear 
span Sn of a finite number of eigenmodes: 

N 

g = k ™<7™ ( 4 ) 
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Figure 1: The part of the spectrum for acoustic modes. The set of neutral eigenmodes outside of the corotation 
interval with u> < mQ(r2) and some growing eigenmodes with u;, > 10 -4 f2o and to r < mflo can be seen, q ~ 1.58, 
w = 300%, 7 = 5/3, Tn — 5. 



The system {q n } forms the basis in this linear subspace and each basis vector has the form: 

q n = K, v v , h}n x = {q\ q 2 , q 3 } n x e mifi , 

where {v r ,v v , h} n are solutions to the one-dimensional boundary problem described in the previous section. We 
assume that the coordinates of q in Sn reflect the temporal behaviour of q: 

where k° are arbitrary complex numbers and uj n are eigen- frequencies of the corresponding eigenmodes q n . 
Thus, the temporal evolution of q is dictated by the matrix expression for its coordinates: 

k = e At K°, (5) 

where k° = (n®, n®, k%) t , k = (ki, kn) T and the exponential is the propagator acting on the initial pertur- 

bation; evidently, A = diag{—iu>i, —iu>2, —iujn}- 

The next fundamental step is to introduce the inner product to measure the specific perturbation magnitude and 
its evolution in time. This is done such that the norm equals to the total acoustic energy of the perturbation: 

Iff Sh 2 \ 
Ea= 2 / 9 V r + * + ) rdrAk P ( 6 ) 

In a bounded flow the latter satisfies the equation [13] : 



— r~ = — f r ~r pSvrSv^r dr dip (7) 
dt dr 



Sound waves cannot carry energy out of the fluid so the integral of the acoustic energy flux across the flow boundary 
vanishes and E a changes only due to the work done by the Reynolds forces. Consequently, it seems natural to define 
the perturbation norm via ©, since the evolution of E a (t) reflects how much energy a given perturbation obtains 
from the basic flow. Next, 

' ' \(Re(q 1 e imv )) 2 + (Re(q 2 e im ' p )) 2 + 



E " 2 



+ -i (Re( q 3 e" n ^)f 
a 1 



rdrdip = 



4 



1-2 



n 



1*1 9*9 1 q * o 

q L q L +q z q z + —q 6 q 6 

or 



rdr, 



where the asterisk denotes the complex conjugation. 

Hence, the inner product between two arbitrary vectors of the linear span (/, g) must have the form: 



(/,s)=-/p[/v+/ 2 y+^/ 3 v 

ri 



rdr (8) 



Finally, the vector norm, i.e. the total perturbation energy is 

N 



{q,q) =E a = \\k\\ 2 = Kt^jMij = K+Mk (9) 



with My = (q~i, q~j) being the metrics in Sn; Here f stands for the Hermitian transpose. 
To conclude this section, we define the growth factor as: 



\\K(t)f ||e A ^°p 

9 ^ = luoW = 11.0112 ( 10 ) 



3.2 The optimal growth 



An instability study looks for an initial perturbation with maximum growth in a finite time interval. This is a kind 
of variational problem when one should fix the time interval T and find the most suitable k° giving the highest value 
of g(T). The mathematical formulation of this problem has been in common use in the hydrodynamical literature 

We need to calculate the propagator norm (see eq. © ), since by definition it is equivalent to the optimal growth 
G (cf. {TOD ): 

G(T)= max 3 (T) ^ ||e AT || 2 (If) 



The matrix M is positive definite, thus the Cholesky decomposition is permissible: M = F^F. This means that we 
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can transit to the Euclidean metrics in Sn- \\k\\ 2 = ||Fk|||. Thus 



e At K°\\ 2 HFe At K°||l 



G(T) = max — - — 7 -- T : — = max 



k°/o || 1 1 2 kVo ||Fk°||5 

-■^'^yB-n^-B (i2, 

The last term in (fT^|) is actually the two-norm of a N x N complex matrix which is equal to the first (the largest) 
singular value, consequently: G(T) — a\ ^Fe AT F -1 ^. 

The singular value decomposition of Fe AT F _1 gives <j\ and the associated right singular vector Vi. The initial 
unit vector n° = F _1 vi corresponds to the perturbation that attains the largest possible energy growth G at the 
moment T. 



Comment on the numerical method 

We have seen that the optimization procedure reduces to a few numerical steps. First, we choose a finite set of 
eigenvalues w n and precisely calculate the radial profiles of eigenmodes {v r , v v , h} n . The metrics M is determined then 
by the numerical integration according to ([5]) . This is the most time consuming task. When we obtain M the standard 
linear algebra methods allow us to get G(t) as well as g(t) for any perturbation with the given initial condition k . 
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Figure 2: The optimal transient growth for a linear combination of different sets of cigcnmodes. From bottom to up: 
dashed curves show 2,4,10,19 growing eigenmodes, solid curves show 2,4,10 neutral eigenmodes (displayed in Fig. 1). 
The bold curve represents the optimal transient growth for a linear combination of both 19 neutral and 19 growing 
eigenmodes and covers the solid curve that corresponds to 19 neutral eigenmodes. The fixed parameters are the same 
as in Fig. 1. 
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Figure 3: The bold curve is the optimal transient growth for a linear combination of the first 5 neutral eigenmodes with 
u> > mil(ri). Other curves denote the energy behaviour of specific initial perturbations depending on time. These are 
the perturbations that attain the largest possible growth at the moments t marked by corresponding numbers. The 
fixed parameters are the same as in Fig. 1. 
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4 Results 



4.1 Combinations of eigenmodes 

First we would like to clarify which combinations of eigenmodes can achieve the most substantial energy growth. We 
take the part of the spectrum displayed in Fig. 1 and calculate G(t) separately for groups of 2,4,10,19 neutral modes 
when the corotation radius is beyond the outer boundary and for groups of 2,4,10,19 growing modes with u>; L > 10~ 4 
and u) r < mQo (see Fig. 2). Obviously, a number of neutral modes is able to exhibit a much greater energy increase 
rather than the growing modes. Moreover, the unification of all 38 modes does not affect G noticeably. This is 
illustrated by the bold curve. One can also see that G-curve is slightly modified when passing from 10 to 19 neutral 
modes. The numerical test shows that this is the case for larger combinations of neutral modes. Furthermore, various 
combinations of the growing modes with uj r > mflo, i.e. with the corotation in the inner part of the flow, were 
examined and did not reveal any substantial growth. 

The fact that slow neutral modes are of the most interest from the standpoint of transient growth is also supported 
by the calculation of G for the fast neutral modes with lu > mSl(ri), i.e. with the corotation radius inside the inner 
boundary of the flow. The result is shown in Fig. 3. Here the curves g(t) are also added to see the temporal behaviour 
of specific initial perturbations. Mathematically, the weak energy growth is explained by small off-diagonal elements 
in M and, consequently, by an almost normal propagator that governs the initial disturbances. 

If the propagator is fully normal G = g = const = 1. In fact, for the unit k. eq. ()10p gives: 

N 

g(t) = \\e At K°\\ 2 = E <V< t 4 e -^ t M„ fc = 

n— l,k— n 
N 

= J2\<\ 2 C 2UjTtM nn (13) 
n=l 

which is unity for neutral modes and is bounded by the largest increment for growing modes. 

Thus we conclude that the optimal growth arises mostly when considering the first 10-15 slow neutral modes with 
lu < mf2(r2). We keep it in mind and present below the results concerning this kind of eigenmodes. 

Besides, one can notice that there is the typical scale of G variation. Though G starts from unity, it never returns 
back (see Fig. 2) and has extrema separated by the characteristic time interval which is evidently determined by the 
shortest difference w$ — uij that is present in Fe^ T F -1 and equals to the step between the eigenvalues in the spectrum. 
Then, the shorter the variation time, the weaker is the transient growth. Note that this dependence is very steep (cf. 
Fig. 2 and Fig. 3). 

4.2 Parametric study of the optimal growth 

In this section we examine the temporal growth by varying free parameters of the problem, i.e. the azimuthal 
wavenumber m, the flow radial size w, and parameter q which defines the rotational profile and pressure distribution 
in the basic flow. The last one is certainly important since it describes the deviation from the pure Keplcrian 
rotation and hypersonic feature of the flow. The research is carried out for the particular case 7 = 5/3 specifying 
the compressibility of the fluid. We restrict ourselves to study the combination of the first 10 slow neutral modes, as 
discussed above. 

Let us first look what happens when we change m, i.e. the characteristic length scale of perturbations. The result 
is shown in Fig. 4. The optimal growth increases with azimuthal wavenumber starting from m = 2. So perturbations 
with smaller length scale have more capacity to gain energy from the background. Note that the timescale of G 
variations weakly depends on m and is nearly 70f2 ~ 1 or 10To, where To is the Keplerian period at tq. Furthermore, 
the first maximum of G occurs at ~ 6To for the above parameters. The exception is the curve corresponding torn = 1. 
It has a longer timescale of variations and attains higher values of G. 

In Fig. 5 and Fig. 6 we present the optimal growth curves for different w and q. The later G attains its maximum 
value, the greater energy growth of perturbations is possible. As we have noticed before, the large scale variation 
of G(t) is defined by the step in the acoustic spectrum which decreases while passing to more extended flows or 
approaching the Keplerian rotation. The last fact also implies the Mach number growth. The explicit dependence on 
the Mach number M defined in the second section of the paper is shown in Fig. 7. Here G max is the first maximum in 
the optimal growth curve and T max is the time at which G = G max . Note that there are several steep changes in the 
overall monotonic growth of G max . This is due to sudden changes in the set of eigenvalues and eigenfunctions we take 
to analyse the optimal growth. Indeed, by convention we take the first 10 neutral modes with the corotation radius 
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Figure 4: The optimal transient growth for a linear combination of the first 10 neutral eigenmodes with lj < mQfa) 
corresponding to different azimuthal wavenumbers. The dashed curve corresponds to m = 1. For solid curves from 
bottom to up: m = 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 30, 50. q = 1.58, w = 300%, 7 = 5/3. 

beyond the outer boundary, so our frequency frame is tu cr — mf2(r2). At the same time the acoustic spectrum reveals 
the continuous shift with respect to Lu cr due to variation of free parameters. Since the spectrum is discrete there are 
moments when some eigenvalue crosses the point uj cr and appears (disappears) in the spectral part that we examine. 

However, in general view G max and T max are pretty proportional to the Mach number, at least in the range of q 
studied numerically. These dependencies can be approximately written as: 

G max ~ 38 + 5.6 (M - 5) 

T max ~ 5.5 + 1.2 (M - 5) 

Note that here T max is given in units of To for conveniency. Let us remind that for a barotropic flow the Mach number 
approximately equals to the ratio of the radial size of the torus to its half-thickness: M ~ 2(r-2 — r{)/h. We should 
mention here as well that [IT] and [TU] studied the boundary layer and the plane Couette problem (and an unbounded 
constant shear) correspondingly. They included both viscosity and compessibility and also found that the capacity to 
transient growth increases with increasing Mach number although this dependence seems to be less strong than in the 
present research. In fact, this may be due to influence of viscosity. 

This result is quite important, since in astrophysics we are dealing with almost Keplerian and thus highly supersonic 
flows. It is remarkable that hypersonic feature of astrophysical flows was indeed the obstacle in the traditional 
eigenvalue analysis since the well-known acoustic instability vanishes while proceeding to high Mach numbers. In spite 
of this fact the specially selected combination of eigenfunctions is able to reach even the larger energy growth for 
almost Keplerian flows with q — > 3/2. 

However, one must still keep in mind that the latter conclusion is done in two-dimensional approach and must be 
checked for group of eigenmodes accurately calculated in three dimensions. Nevertheless, we hope this remains valid 
since the transient growth fenomenon arises from the main feature of the flow considered, i.e. the shear. 

The last two plots are in a sense an addition to the results obtained. Fig. 8 shows the temporal behaviour of 
individual perturbations which is drawn over by the bold curve that represents the optimal growth. The numbers 
denote the time T when the corresponding curves touch the bold one, i.e. when the specific perturbations reach the 
largest possible growth factor G(T). At last we present Fig. 9 to show the long time behaviour of the optimal growth. 
The overall result is that starting from unity G(t) never returns back despite variations with a notable amplitude - 
the feature that we noticed before when commenting on Fig. 2. Here this is confirmed for a longer time interval. The 
consequence is that one can always choose initial perturbations that will reveal a substantial growth at any moment 
in future. Certainly, the necessary condition for this property to exist is the absence of dissipation in our approach. 



Figure 5: The same as in Fig. 4 but for m = 10 and different values of w. For maxima of G (from bottom to up) 
to = 100%, 200%, 250%, 300%, 500%, 1000%. 




Figure 6: The same as in Fig. 5 but for w = 300% and different values of q. For maxima of G (from bottom to up) 
q = 1.58, 1.56, 1.54, 1.52, 1.51, 1.507. 
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Figure 7: G max and T max depending on Mach number as defined in the text. The numbers denote some of the 
corresponding values of q. w = 300%, m = 10, 7 = 5/3. 
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Figure 8: The bold curve is the optimal transient growth displayed for to = 300%, 7 = 5/3, m = 10 and q = 1.507 in 
Fig. 6. The other curves denote the energy behaviour of a specific initial perturbations depending on time. These are 
the perturbations that attain the largest possible growth at the moments of t marked by the corresponding numbers. 




Figure 9: The same as in Fig. 6 but in a longer time interval. The dotted, dashed and solid curves have been calculated 
accordingly for q = 1.58, 1.51, 1.507. 
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5 Summary 



Our study shows the significance of the transient growth concept in application to astrophysical problems. We 
considered here the dynamics of global perturbations in a sense that free boundary conditions which are natural in 
astrophysical practice were involved. Previous studies of a pure hydrodynamical instability in axisymetric flows with 
free boundaries have been carried out by methods of the traditional cigenanalysis, i.e. by searching for growing modes 
with exponential dependence on time. This way was proved to be unable to solve the problem of angular momentum 
transfer in almost Keplerian flows with finite compressibility, since the surface gravity eigenmodes are damping [13j and 
the growing sonic eigenmodes have negligible increments. However, the present two-dimensional analysis demonstrates 
that situation noticeably changes if one turns to the dynamics of a linear combination of eigenmodes. In fact, we studied 
mainly a finite number of slow neutral eigenmodes from acoustic spectrum of the basic flow with lu < mfl(r2). It turns 
out that such a combination exhibits a substantial transient acoustic energy growth g on the dynamical timescale. 
Its magnitude strongly depends on parameters of the problem, i.e. the azimuthal wavenumber of the modes, the 
radial size and the rotation profile of the background flow. The numerical tests produce g oc 10 2 in a few typical 
Keplerian periods. Possibly the most notable fact is that g grows as we approach the Keplerian rotation. This result is 
opposite to the conclusions of eigenanalysis according to which the sonic instability ceases to affect the basic flow just 
because of the small sound speed. If one introduces the Mach number M as the ratio of the shear velocity difference 
between the boundaries r\ and r-i to the maximum sound speed in the flow, then the maximum optimal growth and 
the corresponding time interval are approximately proportional to M. So the results of the present two-dimensional 
analysis indicate that the flat hypersonic tori in the vicinity of gravitating objects can exhibit a large linear transient 
growth of global perturbations which may be related to the problem of angular momentum transfer. However, the 
conclusions we made here certainly must be checked by complete three-dimensional calculations. 
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